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Abstract 

Several explicit examples of multi-particle quasi exactly solvable 'discrete' quantum 
mechanical Hamiltonians are derived by deforming the well-known exactly solvable 
multi-particle Hamiltonians, the Ruijsenaars- Schneider- van Diejen systems. These are 
difference analogues of the quasi exactly solvable multi-particle systems, the quantum 
Inozemtsev systems obtained by deforming the well-known exactly solvable Calogero- 
Sutherland systems. They have a finite number of exactly calculable eigenvalues and 
eigenf unctions. This paper is a multi-particle extension of the recent paper by one of 
the authors on deriving quasi exactly solvable difference equations of single degree of 
freedom. 



1 Introduction 

Recently a recipe to obtain a quasi exactly solvable difference equation from an exactly 
solvable difference equation is developed by one of the present authors []]. In the present 
paper we apply the recipe to obtain multi-particle quasi exactly solvable difference equations. 
A quantum mechanical system is called Quasi Exactly Solvable (QES), if a finite number of 
eigenvalues and the corresponding eigenfunctions can be determined exactly [2]. Since the 



number of exactly solvable states can be chosen as large as wanted, a QES system could be 
used as a good alternative to an exactly solvable system for theoretical as well as practical 
purposes. Many examples of QES systems of single degree of freedom have been known for 
some time, whereas the known examples of those with many degrees of freedom are rather 
limited in their structure [3j 0]. They are all obtained by deforming Calogero- Sutherland 
systems [SI El 13) the exactly solvable Hamiltonian dynamics based on root systems. Among 
them those based on the BC type root systems (rational and trigonometric) and on the A 
type root systems (trigonometric) can contain an arbitrary number of particles. Therefore the 
examples of multi-particle QES are infinite in number. In the present paper we derive several 
examples of multi-particle QES difference equations by deforming Ruijsenaars-Schneider van 
Diejen (RSvD) systems [HJ [9j [lOj [11] , which are the difference equation analogues of Calogero- 
Sutherland systems. To be more precise, we derive multi-particle QES difference equations 
by deforming RSvD systems based on BC (rational and trigonometric) and A (trigonometric) 
type root systems. Thus the examples of multi-particle QES difference equations are now 
infinite in number. The Hamiltonian of RSvD systems has two kinds of interaction terms, 
the 'single particle interaction' part and the 'multi-particle interaction' part. The latter is 
kept intact and the former, the 'single particle interaction' part is deformed according to the 
recipe given in the recent paper pp. 

This paper is organised as follows. In section [2] the Hamiltonian of the rational BC type 
theory is derived and the finite dimensional invariant polynomial space is identified. The 
trigonometric BC type theory is explained in section[3l The Hamiltonian of the trigonometric 
A type theory and the finite dimensional invariant polynomial space are derived in section 
HI The final section is for a short summary and comments. 

2 Rational BC type theory 

The first multi-particle QES Hamiltonian is a simple deformation of the rational Ruijsenaars- 
Schneider-van Diejen (RSvD) [8j [9] system based on the BC n root system. Here n is the 
rank of the root system as well as the degree of freedom with coordinates and conjugate 
momenta: 

def / \ imn def / \ 

X = {X 1 ,X 2 , ■ ■ ■ ,X n ) e K , p = [Pl,P2, ■ ■ ■ ,Pn>- 
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As always, the momentum operator is realised as a differential operator pj = —idj = 
—id/dxj. In 'discrete' quantum mechanical Hamiltonians the momentum operators appear in 
exponentiated forms e ±Pj = e Tldj , in contrast to the ordinary quantum mechanics, in which 
the momentum operators appear as polynomials. Thus their action on the wavefunction is 
a finite shift in the imaginary direction: 



e ±P:i ip(x) = il>(xi, Xj-x,Xj =F i, x j+1 



, ^j + i , ■ ■ ■ , X n j , 



leading to difference Schrodinger equations. The quasi exactly solvable rational Hamiltonian 
has the following general form [Q [TJl [11] : 



H = Y.(VVj ePj \fr? + \Jv? e ~ Pi VVj- v j -V?) + a M, MeN, (2.1) 

3=1 

n 

= Yl {VVj e ~ idj + >fi? ei9j VVj-Vj- v f) + <*M- (2.2) 



Here c^m is a compensation term indexed by a natural number A4, to be specified shortly in 
(12. 5j) and ( 12.7ft . The potential function Vj has the following general form: 

n 

Vj{x) = w{xj) \\ \\ v( Xj + ex k ), (2.3) 

fc=l £=±1 

Mi 

in which the 'multi-particle interaction' part v is not deformed 

v{y)^l-i 9 -, g>0. (2.4) 

y 

Whereas the 'single particle interaction' part w allows two types of deformation for QES, 
corresponding to the linear and quadratic polynomial deformations introduced in section 3.1 
afffl: 

n 

Type I : w(y) = (a 2 + iy)w {y), a M (x) = M^x 2 j, (2.5) 

3=1 

Type II : w{y) = {a x + iy){a 2 + iy)w {y), (2.6) 



, x def 

a M {x) 



M(M-l + J2 a ^ + 2 (n-l)g)J2 x2 j' ( 2 - 7 ) 

a=l j=l 

with a common undeformed wo(y) [T2l [TTJ 

/ \ def Il!-3( a « + W) n /n oN 

"" (!/)= ■ °° >0 - (2 ' 8) 
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As a single particle dynamics, the above wq corresponds to the deformation of the harmonic 
oscillator with the centrifugal potential. The corresponding eigenfunctions are the Wilson 
polynomials [121 03] • 

The main part of the Hamiltonian is factorised [H [121 E] : 



3=1 



Aj 



def 



a, 



it 



V,e-* 9 * - JVfei 



(2.9) 
(2.10) 



exhibiting the hermiticity (self-adjointness) of the Hamiltonian. 

The pseudo ground state wavefunction O is defined as the one annihilated by all the Aj 
operators: 

-Vo = 0, j = l,...,n. (2.11) 

It is given by 

nL=2 r ( a a+*zi) 



Type I : <f> (x t 



Type II: <f>o(x) 



clef 



clef 



n 



T(2iXj) 
Yll=i r {a a + ixj 



n 



r(2zx j ) 



TT TT T ^ + i ( x 3 + ex k)) 

, J- I ' I • ( 

i<j<fc<n£=±i 

n n 



l<j<Kne=±l 



r(i(xj + £x fc )) 

r(g + i(xj + gXfc)) 
r(i(xj + ex k j) 



(2.12) 
(2.13) 



in which an abbreviation |/| *= V //* is used. It is obvious that these 0o are square integrable 
in the principal Weyl chamber of BC n : 



[ (t) {x) 2 d n x < oo, = {x e W 1 | xi > x 2 > ■ ■■ > x n > 0}, 

Jpw 



(2.14) 



and they have no node or singularity in PW. By the similarity transformation in terms of 
the above pseudo ground state wavefunction <fio(x), we obtain a Hamiltonian Ti. leading to a 
finite difference eigenvalue equation with rational potentials fl2.3l) - fl2.8p [H [12] : 

n 

n = 1 o n o 0o = J2{ v ^ e ~ idj ~ x ) + v i( x Y( ei9j - !)) + «m(4 (2.15) 



7i0 = £0, 0(x) = o (x)Pm(x) 



HP^M = £P 



A4- 



(2.16) 



In the undeformed limit, i.e. w = Wo and a^i = 0, the theory is the exactly solvable 
rational BC n RSvD [9] system with a Hamiltonian 7io- The corresponding 0o becomes the 
true ground state wavefunction. The exact solvability means that Tio maps a Weyl-invariant 
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polynomial in {xj} into another of the same degree. The BC n Weyl-invariant polynomials 
in {xj} are simply symmetric (under any permutation j «-> k) polynomials in {x?}. For later 
convenience, let us introduce a monomial symmetric polynomial 

m x({Vj}) =m( Ali .„ 1 A B )(2/i,...,2/ n ) = f y l i---Vni ( 2 - 17 ) 

(2l,.--,in) 

where the summation with respect to . . . ,l n ) is taken over all distinct permutations of 

In the deformed theory, it is straightforward to demonstrate that 7i maps a symmetric 
polynomial in of degree equal or less than A4 into another: 

nV M CV M , dimV M =( M + n ), (2.18) 

V M = Span[m (/li ... A) ({xf}) | < ^ < M, 1 < j < n] . (2.19) 

This establishes the quasi exact solvability. The proof that an invariant polynomial is mapped 
into another goes almost parallel with the undeformed theory. One simply has to verify the 
vanishing of the residues of the simple poles at qj = ±qk, 0, ±i/2. The other step is to show 
that a symmetric polynomial in with degree m (< M. — 1) is mapped to another with 
degree m + 1, whereas a symmetric polynomial of degree M. remains with the same degree. 
This part goes almost the same as in the single particle case shown in PQ, since it is caused 
by the deformation of w, the single particle interaction part. The present examples are 
the difference equation version of the QES theory called rational BC Inozemtsev systems 
discussed in section 6 of Sasaki- Takasaki paper [3] . 

3 Trigonometric BC type theory 

The next example is a QES deformation of the trigonometric BC n RSvD system. Because 
of the periodicity of the trigonometric potential, we introduce a slightly different notation 
for the dynamical variables: 

e = (0 1} e 2 , ... , e n ) er, o < % < vr, Xj = cos e j: Zj = . (3.1) 

The dynamical variables are 6. We denote D 5 = z^. Then q Dj is a g-shift operator, 

q Dj f(z) = /Oi, ■ • • , qzj, Zj+i, 
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with < q < 1. The quasi exactly solvable trigonometric Hamiltonian has the following 
general form [H [121 E] : 



n 

n = E ( ^ j \A? + ^ v 7 ^ - ^ - + 



ax- 



(3.2) 



The compensation term a M is given in (13. 5p . The potential function Vj consists of the 'single 
particle interaction' part w and the 'multi-particle interaction' part v: 



V 3 {z)^w{z 3 )Y\X\v{z 3 zi), 



fc=l £ = ±1 



dd 1 - a y , . def n \ f \ t \ dcf ri Q =2( 1 _ a all) 
v {y) = , W \V) = I 1 - aiy)w {y), w (y) - 



i-y 



[l-y 2 )(l-qy 



2V 



dcf 



(q M - l)q l a 2 ^ n ^ aia 2 a^a i a 5 S ^{z j + z- : ). 



(3.3) 

(3.4) 
(3.5) 



As before, the 'multi-particle interaction' part v is not deformed but the 'single particle 
interaction' part w is multiplicatively deformed by a linear term from wq. As a single particle 
dynamics, the above wq corresponds to the deformation of the Poschl-Teller potential. The 
corresponding eigenf unctions are the Askey- Wilson polynomials [T2l [T3] . 
The main part of the Hamiltonian is factorised [H [T2J [EE] : 



Ai^iU D \hf*-q-^jvX A) 



and the pseudo ground state wavefunction 



D, 



V* q'^ Dj 



<f>o( z ) 



def 



n 



(z^,q)c 



1 ria=l( a a 2 i i ?)oo 1<? - 



n n 



<j<k<n e=±l 



(zjzl ; q)oc 
{a Q ZjZ% ; q) c 



(3.6) 
(3.7) 

(3.8) 
//. Here the 

standard notation (a ; q)^ = f n^Lo(l — ac l n ) anc ^ l/l ^ V77* are used. It is obvious that O 
has no zero or singularity in the principal Weyl alcove of BC n 

pw t = {e em n \ 71 > e 1 + e 2 > e 1 > e 2 > ■ ■ ■ > e n > 0}, (3.9) 

so long as the parameters are restricted to 



is obtained as the common zero-mode of all Aj operators Aj(f)o = 0, j 



1 < a a < 1, a = 0, 1, . . . , 5. 



(3.10) 



Then the square integrability of </>q, f PW §\d n Q < oo is trivial. By the similarity transforma- 
tion in terms of the above pseudo ground state wavefunction 4>o(z), we obtain a Hamiltonian 
7i leading to a finite difference eigenvalue equation with rational potentials in z and z~ x 

n 

W = 1 o H o O = ^(V,(z)(g^ - 1) + ^(^)*(g-^ - 1)) + ^(4 (3.11) 

3=1 

H4> = £4>, <f>{z) = <p Q {z)P M {z) ^ HP M =£P M . (3.12) 

In the undeformed limit, i.e. w = w and = 0, the theory is the exactly solvable 
trigonometric BC n RSvD [9] system with a Hamiltonian 7^ - The corresponding O becomes 
the true ground state wavefunction. The exact solvability means that 7Yo maps a Weyl- 
invariant polynomial in {xj = cosOj = \{zj + zj 1 )} into another of the same degree. The 
BC n Weyl-invariant polynomials in {xj} are simply symmetric (under any permutation 
j <-> k) polynomials in {xj} or in {zj + zj 1 }. The eigenfunctions of Ti,o are the BC type 
Jack polynomials [14] . 

In the deformed theory, it is straightforward to demonstrate that 7i maps a Weyl-invariant 
polynomial in {xj} of degree equal or less than M. into another: 

HV M C Vm, dimV^ = ( M n +n ), (3.13) 

V M = Span[m (/li ... iin) ({^ + zj 1 }) \ < l s < M, 1 < j < n] . (3.14) 

This establishes the quasi exact solvability. The proof that an invariant polynomial is mapped 
into another goes almost parallel with the undeformed theory. One simply has to verify the 
vanishing of the residues of the simple poles at Zj = z^ 1 , ±1, ±g 1//2 , ±g _1//2 . The other step 
is to show that a symmetric polynomial in {xj} with degree m (< A4 — 1) is mapped to 
another with degree m + 1, whereas a symmetric polynomial of degree Ai remains with the 
same degree. This part goes almost the same as in the single particle case shown in PQ, 
since it is caused by the deformation of w, the single particle interaction part. The present 
example is the difference equation version of the QES theory called trigonometric BC type 
Inozemtsev system discussed in section 7 of Sasaki- Takasaki paper [3]. 
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4 Trigonometric A type theory 



The QES deformation of the trigonometric A n _\ RS [8] system goes almost parallel with the 
previous example, or even simpler. For the A-type theory, it is customary to consider A n ^\ 
and to embed all the roots in M. n . This is accompanied by the introduction of one more degree 
of freedom, 9 n and p n . The genuine A n _\ theory corresponds to the relative coordinates and 
their momenta, and the extra degree of freedom is the center of mass coordinate and its 
momentum. The Hamiltonian takes the general form (13.21) with the potential function Vj 

n 

def 



VjW^wiz^Hviz^ 1 ), (4.i; 



k=l 



def 1 - a y i s def , s 

v \y) = > w (y) = 1 - °>iy, ( 4 - 2 ) 

1-2/ 

n 

<y M (z) = (q M - lK^ai + Z J 1 )- ( 4 - 3 ) 

3=1 

The undeformed theory, the trigonometric A n _i RS system |SJ has w = 1 and = 0. The 

main part of the Hamiltonian is factorised as in (13.61) and (13.71) . The pseudo ground state 
wavefunction 0o annihilated by all Aj reads 

1 TT ( z i z k 1 ; ?)oo 



<f>o( z ) = 



(4.4) 



f = \ {aizj ; 9)00 1 < J 1 < t<„ (ao% 1 ; 

which has no node or singularity in the principal Weyl alcove of A n _\. 

pw t = {e e R n 1 7T > e l > e 2 > ■ ■ ■ > e n > o}, (4.5) 



so long as the parameters are restricted to — 1 < ao, a\ < 1. The similarity transformed 
Hamiltonian Ti, in terms of 0o nas the same form as (13.111) . (I3.12p . 

In the undeformed limit, i.e. w = 1 and = 0, the theory is the exactly solvable 
trigonometric A n _i RS system with a Hamiltonian 7f - The eigenf unctions of 7Yo are the 
well-known Jack polynomials in {zj} [14j . Since all the coefficients of the eigenvalue equation 
HoV 9 = ^-V 9 are r eal, the Jack polynomials in are also eigenf unctions. In other words 

Tio maps a symmetric polynomial in {zj} into another of the same degree. 

The eigenfunctions of the deformed Hamiltonian Ti, are still symmetric polynomials in 
{zj} and {zj 1 }, but truncated to the maximal power of A4 for each variable, due to the 
single particle interaction term w ( 14. 2 p and the compensation term (14.31) . This form of the 
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compensation term is necessary for the hermiticity of the Hamiltonian. In other words 7i 
has the invariant subspace 

HV M <ZV Ml dimV M = ( 2M + n ), (4.6) 

Vm = Span[m 0li ... A) ({ % }) | -M < l 4 < M, 1 < j < n] . (4.7) 

This establishes the quasi exact solvability. The proof that a symmetric polynomial is 
mapped into another goes almost parallel with the undeformed theory. One simply has 
to verify the vanishing of the residues of the simple poles at Zj = Zk- The other step is to 
show that a symmetric polynomial in {zj} with degree m (-A4 + 1 < m < A4 — 1) is mapped 
to another with degree mil, whereas a symmetric polynomial of degree ±Ai remains with 
the same degree. This part goes almost the same as in the single particle case shown in [TJ, 
since it is caused by the deformation of w, the single particle interaction part. The present 
example is the difference equation version of the QES theory called trigonometric A type 
Inozemtsev system discussed in section 8 of Sasaki- Takasaki paper [3j. 

5 Summary and Comments 

Quasi exactly solvable multi-particle difference equations are derived by deforming rational 
and trigonometric BC type RSvD systems as well as the trigonometric A type RS systems. 
The method of multi-particle deformations is a simple extension of the single particle case 
recently developed by one of the authors [TJ. These examples are the difference equation 
analogues of the quasi exactly solvable multi-particle quantum mechanical systems derived 
by Sasaki- Takasaki [3J. 

A few comments are in order. As in the single particle cases, the ranges of the parameters 
can be loosened without losing QES. For example, the six positive parameters in Type II in 
section 2 (12. 6p . a±, . . . , could be replaced by three complex conjugate pairs with positive 
real parts. Among the five real parameters ai,. . . , in (I3.10p . four could be replaced by 
two complex conjugate pairs bi, b{, 62, b\ with modulus less than unity, |&i| < 1, I&2I < 1- 

Let us elaborate on the connection with the trigonometric Inozemtsev systems mentioned 
at the end of section [3J and HI In fact the examples in section [3J and |4] reduce to the Inozemtsev 
systems given in [3J in a certain limit (q — > 1). Let us introduce a parameter c, and set q = 

e -2/C ; aQ = q 9 7 a l = l- Q = 2 ^new_ Moreover get ftnew = ^^^^^-1^-1/2 a ~{n-l) 1^ 
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and (a 2 , a 3 , a 4 , a 5 ) = (q 91 , q 92+1/2 , -q 9 '\ -q 9 '^ 1 ' 2 ) for the BC type theory, ft new = a (n ~ 1)/2 H 
for the A type theory (see [EE]). Then in the c — > oo limit, the Hamiltonian |c 2 7i new reduces 
to the trigonometric Inozemtsev systems discussed in section 7 and 8 of [3]. 

The corresponding statement for the rational theory is somehow complicated and it 
requires a double limit. For the parameters in section [2J let us set (a^^a^a^^ao) = (~v/^p + 
±^{^ _ i) c /^r+ _L^(^ _ l) m o 2 + I) and x j = cxf w . Moreover set H new = 
4(a2a-3a4)~ 1 7i and a 2 = c 2 for the type I theory, 7i new = 4(aia2a3Ci4) _1 7i and a\ = a 2 = 
2c 2 /ui for the type II theory (see [H]). Then a double limit lim ( lim ^c 2 7^ new ) gives the 
Hamiltonian of the rational Inozemtsev system discussed in section 6 of [5J. 

It is interesting to note that the weight function <Pq(x) for the polynomial eigenfunctions 
{Pm(x)} is the zero mode (stationary distribution) of the corresponding deformed Fokker- 
Planck equation |15j . 
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